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Abstract 

We describe the isolation, biological and genetic characterization of a host-range variant of bovine coronavirus (BCoV) detected in water 
buffalo (Bubalus bubalis). By conventional and real-time RT-PCR assays, the virus was demonstrated in the intestinal contents of two 20-day-old 
buffalo calves dead of a severe form of enteritis and in the feces of additional 17 buffalo calves with diarrhea. Virus isolation, hemagglutination 
and receptor-destroying enzyme activity showed that the buffalo coronavirus (BuCoV) is closely related to BCoV but possesses some different 
biological properties. Sequence and phylogenetic analyses of the 3' end (9.6 kb) of the BuCoV RNA revealed a genomic organization typical of 
group 2 coronaviruses. Moreover, the genetic distance between BuCoV and BCoV was proven to be the same or even higher than the distance 
between other ruminant coronaviruses and BCoV. In conclusion, our data support the existence of a host-range variant of BCoV associated with 
enteritis in buffaloes. 

© 2007 Elsevier Inc. All rights reserved. 

Keywords: Buffalo; Enteritis; Bubaline coronavims; Biological characterization; Molecular characterization 


Introduction 

Coronaviruses (CoVs) (order Nidovirales , family Coronavir- 
idae) are enveloped, positive-sense single-stranded RNA particles 
that are responsible for enteric and/or respiratory disease in 
mammals and birds. Their ability to evolve through genetic 
recombination and/or point mutation is well-known so that they 
give rise to new viral genotypes or mutants with different tissue or 
host tropism. Currently, CoVs are organized into three antigenic 
groups with group 2 including bovine-like (subgroup 2a) and 
SARS-like (subgroup 2b) viruses. Bovine coronavirus (BCoV) 
belongs to subgroup 2a together with mouse hepatitis virus, 
sialodacryadenitis virus, porcine hemagglutinating encephalo¬ 
myelitis virus, human coronavirus (HCoV) OC43, human enteric 
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coronavirus (HECV) 4408 (Enjuanes et al., 2000) and the newly 
recognized equine coronavirus (ECoV) (Guy et al., 2000), HCoV- 
HKU1 (Woo et al., 2005), and canine respiratory coronavirus 
(CRCoV) (Decaro et al., 2007a; Erles et al., 2003). There are 
multiple genetic and antigenic evidences that several subgroup 2a 
CoVs, such as HCoV-OC43, HECV-4408, PHEV and CRCoV, 
have arisen as a consequence of trans-species infections caused by 
BCoV (Zhang et al., 1994; Vijgen et al., 2005, 2006; Erles et al., 
2007). BCoV can cause severe diarrhea in newborn calves 
(Snodgrass et al., 1986), “winter dysentery” in adult cows (Cho 
et al., 2000; Saif et al., 1991) and respiratory tract illness in calves 
and cows (Lathrop et al., 2000; Storz et al., 2000). The same virus 
strain could be responsible for simultaneous appearance of enteric 
and respiratory disease in the same animals (Chouljenko et al., 
2001) as well as in both calves and cows (Traven et al., 2001). 
BCoV infection is responsible for severe economic losses in cattle 
herds due to decreased milk production in dairy cows and reduced 
weight gain in beef calves (Saif et al., 1998), with generally low 
mortality rates. 
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Water buffalo ( Bubalus bubalis) represents a ruminant spe¬ 
cies important in the economy of several countries, including 
Brazil, India, Vietnam and some regions of southern Italy, 
where the breed Mediterranean Italian buffalo produces high 
quality milk employed for production of the buffalo “mozza¬ 
rella cheese” (Romano et al., 2001). Recently, bovine-like CoVs 
have been identified in wild or domesticated ruminants, 
including several species of deer, waterbuck antelope (Tsune- 
mitsu et al., 1995), giraffe ( Giraffa Camelopardalis) (Hasoksuz 
et al., 2007), alpaca ( Lama pacos) (Jin et al., 2007) and sable 
antelope ( Hippotragus niger) (Spiro D., R. Halpin, S. Wang, M. 
Hasoksuz, X. Zhang, K. Alekseev, A. Vlasova, D. Janies, E. 
Ghedin, and L. Saif, unpublished). However, Co Vs have never 
been isolated from buffaloes, although there is a single report 
on the detection of BCoV antibodies in Bulgarian buffaloes 
by hemagglutination inhibition and virus neutralization tests 
(Muniiappa et al., 1985). In the same report, hemagglutinating 
activity which could be suppressed by a BCoV-specific bovine 
serum was also demonstrated in some bubaline fecal samples. In 
this paper, we describe the isolation and biological and genomic 
characterization of a bovine-like coronavirus from a buffalo 
herd affected by severe diarrhea and calf mortality. 

Results 

Clinical outbreak 

The outbreak occurred between October 2006 and April 
2007 in a herd of Mediterranean Italian buffalo ( B . bubalis ) in 
Campania (southern Italy). No cattle, sheep or goat farms were 
present in the vicinity of the outbreak. At that time the herd 
consisted of 460 buffaloes, including 215 lactating cows, all 
vaccinated against colibacillosis, clostridiosis and salmonello¬ 
sis. Buffalo calves were removed from their dams shortly after 

Table 1 


Bovine-like CoV RNA titers detected by real-time RT-PCR in buffalo calves 
with diarrhea 


Prot. no. 

Age of the 
animals 

Sample 

type 

RNA copies/pl of 
template 

179/07-A 

20 days 

Intestinal content 

2.13x 10 3 

179/07-B 

20 days 

Intestinal content 

9.41 x 10 2 

179/07-1 

8 days 

Feces 

1.25 x 10 3 

179/07-2 

2 months 

Feces 

4.10x 10 4 

179/07-3 

2 months 

Feces 

5.98 x 10 6 

179/07-4 

5 days 

Feces 

4.55 x 10 4 

179/07-5 

1 month 

Feces 

8.14x 10 4 

179/07-6 

2 months 

Feces 

3.27 x 10 5 

179/07-7 

5 days 

Feces 

2.90 x 10 2 

179/07-8 

12 days 

Feces 

6.71 x 10 3 

179/07-9 

2 months 

Feces 

2.76 x 10 3 

179/07-10 

2 months 

Feces 

7.18x 10 2 

179/07-11 

3 months 

Feces 

5.23 x 10 7 

179/07-12 

1 month 

Feces 

2.04 x 10 2 

179/07-13 

1 month 

Feces 

1.77xl0 5 

179/07-14 

2 months 

Feces 

3.34 x 10 4 

179/07-15 

1.5 months 

Feces 

4.29 x 10 2 

179/07-16 

2 months 

Feces 

6.83 x 10 3 

179/07-17 

2 months 

Feces 

2.69 x 10 4 


A 



Fig. 1. Virus isolation of the bubaline CoV on human rectal tumor cells. 
(A) Cytopathic effect (syncytia) caused by CoV strain 179/07-11 (72 h 
postinoculation). (B) Cytoplasmic fluorescence detected by the immunofluores¬ 
cence assay using a BCoV-specific serum (24 h postinoculation). 

their birth and placed in separate hutches according to the 
gender where they were hand-fed fresh colustrum for 5 days. 
Neonatal mortality was firstly observed in October 2006 in 30 
5- to 20-day-old calves (out of 40 newborns) that displayed 
severe diarrhea and died despite treatment with antibiotics 
(oxytetracyclin and amoxicillin). Simultaneously, gastroenteric 
disease was also observed in older calves (1-3 months of age). 
On April 2007, neonatal mortality and enteric signs in calves 
persisted and two carcasses of 20-day-old dead calves together 
with fecal samples from additional 17 diseased calves were 
submitted to our laboratory for routine analysis. 

Identification of a bovine-like CoV in buffalo calves 

At necropsy, the carcasses of the two dead buffalo calves 
showed severe gastroenteritis, with enlargement of the mesen¬ 
teric lymph nodes and gall bladder. By conventional RT-PCR 
(Erles et al., 2003), bovine-like CoV RNA was detected in the 
intestinal content of the dead animals as well as in all 17 fecal 
samples from calves with diarrhea. Sequence analysis of the S- 
gene fragments showed a 100% nucleotide identity among the 
CoV strains detected in different animals. 
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Table 2 


Hemagglutination (HA) and receptor-destroying enzyme (RDE) titers of bubaline CoV 179/07-11 and BCoV strains 


CoV strain 

Origin 

HA titer 3 

4 °C 


37 °C 


RDE titer 3 


Mouse 

Chicken 

Mouse 

Chicken 

Mouse 

Chicken 

BuCoV-179/07-11 

Feces 

512 

— 

128 

— 

128 

— 

BCoV-339/06 

Feces 

128 

64 

64 

32 

64 

16 

BCoV-43 8/06-2 

Nasal swab 

512 

256 

128 

16 

128 

64 

BCoV-9WBL7 

Unknown 

1024 

256 

1024 

64 

256 

— 


a HA and RDE titers per 50 pi of fecal suspension are expressed as reciprocals of the highest dilutions producing HA and resulting in complete disappearance of HA, 
respectively. Titers <2 are indicated as negative (-). 


Gastrointestinal parasites of ruminants, including Cryptos¬ 
poridium spp., were not detected by zinc sulfate flotation or 
Ziehl Nielsen staining, whereas by bacteriological investigations 
only common bacteria with poor pathogenic significance 
(E. coli , Streptococcus spp., Staphylococcus spp., Corynebac- 
terium spp.) were isolated from the intestinal contents or fecal 
samples. Using real-time RT-PCR, bovine-like CoV RNA was 
detected at low titers in the intestinal contents of the dead 
animals, whereas higher viral loads were found in the fecal 
samples of calves with diarrhea, with a peak of 5.23 x 10 RNA 
copies/pl of template in calf 179/07-11 (Table 1). 

Biological characterization of the bubaline CoV 

The fecal sample of calf 179/07-11 was used for virus isolation 
attempts on human rectal tumor (HRT-18) and Madin Darby 
bovine kidney (MDBK) cells. Cytopathic effect (CPE) consisting 
of syncytia formation and subsequent cell lysis was evident at the 
2nd passage on HRT-18 cells (Fig. 1A), where the growth of a 
bovine-like CoV was confirmed by the cytoplasmic fluorescence 
detected by the IF assay using a BCoV-specific serum (Fig. IB). 
At the 5th passage on HRT cells the viral titer reached values of 
10 6 ' 25 TCID 5O /50 pi of viral suspension and 1.97 x 10 8 RNA 
copies/ml of template, as calculated by the Karber’s method and 
real-time RT-PCR, respectively. Replication of strain 179/07-11 
in MDBK cells was not as efficient as in HRT-18 cells, 
considering that CPE was not observed on MDBK monolayers 
and only few positive cells were detected by the IF assay at the 1 st 
passage, with viral titers (as calculated by real-time RT-PCR) 
decreasing progressively in the subsequent passages until defi¬ 


nitive loss of the replication ability. MDBK cells were confirmed 
to be poorly permissive to 179/07-11 replication even when the 
HRT-18 cell-cultured virus was inoculated. 

The isolated strain was evaluated by assessment of hemag¬ 
glutination (HA) and receptor-destroying enzyme (RDE) 
activities in comparison with the BCoV isolates 339/06 (Decaro 
et al., 2007b), 438/06-2 (Decaro et al., in press) and 9WBF7 
(kindly provided by Dr Paolo Cordioli, Istituto Zooprofilattico 
Sperimentale della Fombardia e dell’Emilia Romagna, Brescia, 
Italy). All isolates agglutinated mouse erythrocytes at titers 
higher at 4 °C than at 37 °C. The bovine isolates caused HA with 
chicken erythrocytes at low titers at both 4 °C and 37 °C, 
whereas the bubaline strain 179/07-11 did not hemagglutinate 
those cells despite the high HA titers obtained with mouse red 
blood cells (Table 2). Using mouse erythrocytes, RDE activity 
(loss of the HA pattern formed at 4 °C) was observed for the 
bovine and bubaline isolates (Table 2). 

Genetic characterization of the bubaline CoV 

A 9.6-kb region encompassing the entire 3' end of the viral 
genome (from the 32-kDa to nucleocapsid protein genes) was 
determined through PCR amplification and subsequent sequenc¬ 
ing of overlapping fragments. At the 3' end of viral RNA, 
bubaline CoV strain 179/07-11 had the same genomic organiza¬ 
tion of other ruminant CoVs (Fig. 2). Nine ORFs were identified 
by both ORF Finder program and sequence comparison with 
reference ruminant CoV sequences. All predicted ORFs but the 
4.8-kDa gene were preceded by a repeated intergenic sequence, 
CUAAAC or CCAAAC, which is predicted to interact with the 
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84 
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Fig. 2. Schematic comparison of the genomes (3' end) of different ruminant CoVs. RdRp, RNA dependent-RNA polymerase. Below the diagram, the length in amino 
acids is reported for the encoded proteins of bubaline coronavirus 179/07-11 (BuCoV), bovine coronavirus Mebus (BCoV), giraffe coronavirus US/OH3/2003 
(GiCoV), alpaca coronavirus (ACoV) and sable antelope coronavirus US/OH1/2003 (SACoV). 
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Table 3 


Coding potential and putative transcription regulatory sequences of the 3' end of 
the BuCoV genome 


Putative gene 


Putative TRS 


Gene segment 

Coding sequence 

Start nt position 

TRS sequence 

32 kDa 

200-1036 

187 

CUAAAC 

HE 

1048-2322 

1033 

CUAAAC 

S 

2337-6428 

2330 

CUAAAC 

4.9 kDa 

6418-6549 

6095 

CCAAAC 

4.8 kDa 

6585-6722 

— 

Not detected 

12.7 kDa 

6806-7135 

6725 

CCAAAC 

E 

7122-7376 

6993 

CCAAAC 

M 

7390-8082 

7381 

CCAAAC 

N 

8092-9438 

8079 

CUAAAC 

I 

8153-8776 

8079 

CUAAAC 


viral polymerase during the discontinuous transcription of the 
negative strand subgenomic RNA of the Nidovirales members 
(Table 3). 

Table 4 shows the unique aa changes encountered in the 
encoded proteins of strain 179/07-11 with respect to other ru¬ 
minant Co Vs. The predicted 32-kDa nonstmctural protein (nsp) 
had the same length (278 aa) as most other bovine-like CoVs, 
including ruminant CoVs giraffe coronavirus (GiCoV), alpaca 
coronavirus (ACoV) and sable antelope coronavirus (SACoV), as 


well as CRCoV and HCoV-OC43. Three aa changes were found 
to be unique to strain 179/07-11. The highest aa identity (98.9%) 
of the 32-kDa protein of strain 179/07-11 was found to enteric and 
respiratory BCoVs, as well as to ACoV. The hemagglutinin- 
esterase (HE) protein (424 aa) showed three unique aa mutations 
in comparison with other ruminant CoVs and it was more closely 
related (99% aa identity) to BCoVs DB2 and ENT and SACoV. 
By analysis with the NetNglyc server, the HE of strain 179/07-11 
was found to contain nine potential N-glycosylation sites, ana¬ 
logously to BCoV reference strains. Moreover, the predicted site 
for neuraminidate-O-acetylesterase activity, FGDS, was detected 
at the N-terminus. The spike (S) protein was long 1363 aa and 
matched the best aa identity (98.1%) to BCoV-E-AH65 and 
SACoV-US/OHl/03. Twenty-one potential N-glycosylation 
sites were identified throughout the protein, one more than other 
ruminant CoVs, due to the presence of an additional N- 
glycosylation site, NLS, at position 1260. No deletions or insertions 
were observed in strain 179/07-11 compared with bovine-like 
CoVs. Also, the aa stretch KRRSRR, responsible for the proteolytic 
cleavage of the S protein at residue 768 into subunits SI and S2, 
was conserved. In comparison to mminant CoV reference strains, 
16 unique residues were identified in the S protein of strain 179/07- 
11. Nsp 4.9 kDa was 43-aa long as for BCoVs Mebus, Quebec and 
DB2 and HECV-4408, whereas in other mminant CoVs it is 


Table 4 


Unique amino acid residues observed in structural and nonstmctural proteins encoded by the 3' end of the bubaline CoV genome in comparison with other mminant 
CoVs 


CoV strain 

32 kDa a 



HE 


S 


















205 

217 

275 

108 158 210 

143 

147 

151 

157 

257 

318 

328 

525 

608 

785 

805 

828 

888 

909 

927 

1260 

BuCoV-179/07-11 

K 

N 

F 


H 

F 

> N 

Y 

F 

I 

V 

N 

T 

F 

Y 

G 

K 

D 


T 

K 

R 

A 

N 

BCoV-Mebus 

M 

K 

V 


Y 

8 

; k 

H 

L 

L 

F 

T 

F 

L 

H 

D 

N 

E 


A 

N 

K 

S 

D 

BCoV-ENT 

M 

K 

V 


Y 

5 

; k 

H 

L 

L 

I 

T 

F 

L 

H 

D 

N 

E 


A 

N 

K 

S 

D 

BCoV-Quebec 

M 

K 

V 


Y 

8 

; k 

H 

L 

L 

F 

T 

F 

L 

H 

D 

N 

E 


A 

N 

K 

s 

D 

BCoV-DB2 

M 

K 

V 


Y 

5 

; k 

H 

L 

L 

I 

T 

I 

L 

H 

D 

N 

E 


A 

N 

K 

s 

H 

GiCoV-US/OH3/2003 

M 

K 

V 


Y 

8 

; k 

H 

L 

L 

F 

T 

F 

L 

H 

D 

N 

E 


A 

N 

K 

s 

D 

ACoV 

M 

K 

V 


Y 

5 

; k 

H 

L 

L 

I 

T 

I 

L 

H 

D 

N 

E 


A 

N 

K 

s 

D 

SACoV-US/OH 1/2003 

M 

K 

V 


Y 

8 

; k 

H 

L 

L 

F 

T 

F 

L 

H 

D 

N 

E 


A 

N 

K 

s 

D 

CRCoV-4182 

M 

K 

V 


Y 

5 

; k 

H 

L 

L 

I 

T 

I 

L 

Y 

G 

N 

E 
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s 

D 

HECV-4408 

NA 

NA 

NA 

Y 

8 

; k 

H 

L 

L 

F 

S 

F 

L 
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G 

N 

E 
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N 

R 

A 

D 

HCoV-OC43 

M 

K 

V 


Y 

A R 

R 

S 

V 

V 

T 

I 

L 

D 

D 

N 

E 


A 

N 

K 

s 

D 

PHEV-VW572 

M 

NP 

NP 

Y 

F 

> K 

H 

F 

L 

F 

S 

V 

L 

T 

D 

N 

E 


T 

N 

R 

s 

D 

CoV strain 

4.9 kDa b 




4.8 kDa b 






12.7 

kDa 

M 


N 




F 





24 

32 


40 


2 

9 

16 

28 

29 

40 

41 

70 


58 

158 

152 

354 

56 

116 

122 

132 

BuCoV-179/07-11 

P 

A 


T 


Q 

D 

K 

G 

F 

N 

W 

F 


L 

S 

H 


S 


K 

I 

S 

M 

BCoV-Mebus 

L 

E 


M 


p 

G 

M 

C 

S 

T 

c 

S 


F 

G 

Q 


N 


R 

T 

L 

R 

BCoV-ENT 

L 

NP 


NP 


p 

G 

M 

V 

S 

T 

c 

S 


F 

G 

Q 


N 


R 

T 

L 

R 

BCoV-Quebec 

L 

E 


M 


NP 

G 

M 

c 

s 

T 

c 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

BCoV-DB2 

L 

E 


M 


P 

G 

M 

V 

s 

T 

c 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

GiCoV-US/OH3/2003 

L 

NP 


NP 


P 

V 

M 

V 

s 

NP 

NP 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

ACoV 

L 

NP 


NP 


P 

V 

M 

V 

s 

T 

s 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

SACoV-US/OH 1/2003 

L 

NP 


NP 


P 

V 

M 

V 

s 

NP 

NP 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

CRCoV-4182 

NP 

NP 


NP 


NP 

NP 

NP 

NP 

NP 

NP 

NP 

s 


F 

G 

Q 


N 


R 

T 

L 

R 

HECV-4408 

L 

E 


T 


Q 

D 

M 

G 

s 

N 

W 

s 


L 

G 

Q 


N 


R 

T 

L 

R 

HCoV-OC43 

NP 

NP 


NP 


NP 

NP 

NP 

NP 

NP 

NP 

NP 

s 


F 

G 

Q 


N 


NP 

T 

L 

R 

PHEV-VW572 

F 

NP 


NP 


P 

NP 

NP 

NP 

NP 

NP 

NP 

s 


F 

G 

Q 


N 


R 

T 

L 

R 


NA, not available; NP, not present. 

a Nsp 32 kDa of HECV-4408 is not available in the GenBank database, whereas it is tmncated in PHEV-VW572. 
b Nsp 4.9 kDa and 4.8 kDa are not present in the CRCoV-4182 and HCoV-OC43 genomes. Nsp 4.9 kDa is truncated in PHEV-VW572. 
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truncated in the C-terminus (only 29 aa in length). The length of nsp 
4.8 kDa was 45 aa as for HECV-4408 and ruminant Co Vs with the 
exception of GiCoVand SACoV (38 aa). Nsp 12.7 kDa showed the 
same size as for HECV-4408 and mminant Co Vs (109 aa), 
although BCoV Quebec exhibits a tmncate form of the protein (87 
aa). In the three nonstructural proteins, strain 179/07-11 displayed 
the highest identity to HECV-4408, with three, seven and one 
unique residues in nsp 4.9 kDa, 4.8 kDa and 12.7 kDa, respectively. 
The envelope (E) protein was 84-aa long exactly as BCoV and 
GiCoV but 1-aa shorter and 2-aa longer than the analogous protein 
of ACoV and SACoV, respectively. No residue unique to the 
bubaline strain was detected at this level. A 100% aa identity in the 


E protein was found between the bubaline strain and GiCoV and 
most BCoV strains, with the exception of BCoV-Mebus (98.8%), 
whereas the bubaline isolate was slightly less related to ACoV and 
SACoV. The membrane (M) protein was 220-aa long and contained 
two unique residues, with the 15 8L change shared with HECV- 
4408. In this protein the highest identity (99.5%) was matched to 
HECV-4408. One potential N-glycosylation and four potential O- 
glycosylation sites were detected at the N-terminus of the M 
protein, as for all bovine-like Co Vs including the closely related 
HECV-4408, with which the bubaline isolates shared an aa 
mutation with respect to other ruminant Co Vs. However, both 
programs indicated that M gene sequence of the ruminant Co Vs 
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Fig. 3. Neighbor-joining trees based on the hemagglutinin-esterase (A) and membrane (B) proteins of group 2 CoVs. For phylogenetic tree construction, the following 
CoV strains were used (GenBank accession numbers are reported in parentheses): BCoV-Mebus (U00735), Quebec (AF220295), DB2 (DQ811784), ENT 
(AF391541), LUN (AF391542), E-AH65 (EF424615), R-AH65 (EF424617), E-AH65-TC (EF424616), R-AH65-TC (EF424618), E-AH187 (EF424619), R-AH187 
(EF424620); GiCoV-US/OH3/2003 (EF424623); ACoV (DQ915164); SACoV-US/OH 1/2003 (EF424621); CRCoV-4182 (DQ682406); HCoV-OC43 (NC_005147); 
HECV-4408 (HE, L07747; M, AY316299); PHEV-VW572 (DQ011855); MHV-A59 (AY700211); SDAV (AF207551); HCoV-HKUl (NC_006577). A statistical 
support was provided by bootstrapping over 1000 replicates. The scale bars indicate the estimated numbers of amino acid substitutions per site. 
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might not contain a signal peptide. Proteins without signal pep¬ 
tides are unlikely to be exposed to the O-glycosylation and N- 
glycosylation machineries and thus may not be glycosylated in vivo 
even though they contain potential motifs. The nucleocapsid (N) 
protein of strain 179/07-11 had a length of 448 aa and was closely 
related (99.3% of aa identity) to BCoV-DB2. Two unique aa 
changes were found in comparison to other bovine-like viruses. The 
I protein (207 aa) encoded by the N-gene internal ORF showed four 
unique residues and the best aa identity (96.1%) to BCoV reference 
strains Mebus, Quebec and DB2. When the 9.6-kb sequence 
corresponding to the 3' end of CoV genome was analyzed, the 
bubaline strain showed over 97.6% nt identity to BCoV strains and 
other mminant Co Vs, whereas it was less genetically related to 
other bovine-like CoVs, including CRCoV (96.5% nt identity), 
HCoV-OC43 (91.4%), HECV-4408 (85.1%) and PHEV-VW572 
(83.1%). 

Phylogenetic analysis with the structural protein HE showed 
that the bubaline strain 179/07-11 clusters together with the 
bovine-like Co Vs, being more related to ruminant viruses 
(Fig. 3). Similar phylogenetic trees were generated for the other 
proteins (data not shown) but the M protein, where strain 179/ 
07-11 forms a separate cluster with HECV-4408 into the bovine 
subgroup (Fig. 3). 

Evaluation of the genetic distance among ruminant CoVs 

In order to better define the taxonomic position of the bubaline 
CoV 179/07-11, a sequence comparison between BCoV reference 
strain ENT and ruminant CoVs GiCoV, ACoV, SACoV and the 
bubaline strain was conduced in the nonstmctural and structural 
proteins (Table 5). For such a comparison, strain ENT was chosen 
among the available BCoV sequences as it is an enteric strain 
isolated very recently (Chouljenko et al., 2001). The results 
showed that strain 179/07-11 is less related to BCoV-ENT than to 
other ruminant CoVs in all major structural and nonstmctural 
proteins, whereas closer relationship (98.8% aa identity) to 
BCoV-ENT with respect to SACoV (97.6% aa identity) was 
found only in the E protein. Analysis of the nt sequence cor- 

Table 5 


Sequence identity (%) of ruminant CoVs to reference BCoV-ENT in 
nonstmctural and structural proteins and in the full-length sequence of the 
genomic RNA 3' end (9.6 kb) 


Gene 

segment 

CoV strain 




BuCoV 

179/07-11 

GiCoV 

US/OH3/2003 

ACoV 

SACoV 

US/OH1/2003 

32 kDa a 

98.2 

98.5 

99.2 

98.5 

HE a 

99.0 

100 

99.7 

100 

S a 

98.0 

98.8 

99.1 

99.3 

4.9 kDa a 

55.8 

100 

100 

100 

4.8 kDa a 

75.5 

82.2 

93.3 

82.2 

12.7 kDa a 

98.1 

99.0 

99.0 

98.1 

E a 

98.8 

100 

98.8 

97.6 

M a 

98.2 

99.5 

99.5 

99.5 

N a 

98.4 

99.5 

99.7 

99.5 

r 

95.1 

98.0 

99.0 

98.5 

9.6-kb 2' end b 

97.8 

99.1 

99.3 

99.2 


a Amino acid identity. 
b Nucleotide identity. 


responding to the entire 3' end of CoV genome confirmed that the 
bubaline strain has a lower genetic relationship (97.8% nt identity) 
to the reference BCoV than other recently isolated mminant CoVs 
(99.1-99.3% nt identity). Confirmation of the conclusions drawn 
from Table 5 was derived with the use of SimPlot analysis. 

Discussion 

In the present report, a bovine-like CoV in water buffalo ( B . 
bubalis) has been described, along with its biological and 
genetic characterization. The bubaline strain was detected in the 
feces of buffalo calves with severe gastroenteritis, causing the 
death of at least 30 animals. A vims (strain 179/07-11) anti- 
genically related to BCoV was isolated on HRT-18 cells, 
producing CPE typical of BCoV and being recognized by a 
BCoV-specific semm in an IF assay. Analysis of the biological 
properties of the bubaline strain showed that, unlike classical 
BCoV strains (Fukutomi et al., 1999), the vims was not able to 
replicate efficiently on MDBK cells or cause hemagglutination 
with chicken erythrocytes. However, Benfield and Saif (1990) 
demonstrated that winter dysentery strains of BCoV from adult 
cattle were also unable to replicate in MDBK cells. Among wild 
mminant CoVs, only few strains have been characterized at a 
biological level. GiCoV displayed HA activity also with chick¬ 
en erythrocytes at a temperature of 4 °C, whereas no data 
are available about its replication efficiency on MDBK cells 
(Hasoksuz et al., 2007). Using chicken erythrocytes, other 
mminant CoVs, including the waterbuck, sambar deer and 
white-tailed deer isolates and some BCoV strains, exhibited 
very low HA titers at 4 °C and no HA activity at 37 °C. Thus, 
the loss of ability to replicate on MDBK cells and to agglutinate 
chicken erythrocytes appears as characteristic of the bubaline 
CoV with respect to the other mminant CoVs. 

Sequence analysis of the 3 ' end of the viral genome showed 
that strain 179/07-11 has a genomic organization similar to 
BCoV, including the presence of three ORFs encoding for small 
nonstmctural proteins between the S and E genes. Moreover, the 
vims has a close genetic relatedness to the bovine CoV sub¬ 
group in the major stmctural and nonstmctural proteins. How¬ 
ever, as shown by sequence and SimPlot analyses, the genetic 
distance between the bubaline vims and reference BCoV strain 
ENT was approximately the same or even higher than observed 
between the BCoV strain and other mminant CoVs. Unique aa 
substitutions were detected in all proteins of the bubaline vims 
but they were mainly within the S protein (Table 4). As ob¬ 
served in BCoV (Chouljenko et al., 1998), GiCoV (Hasoksuz 
et al., 2007) and ACoV (Jin et al., 2007), these changes occurred 
mostly in the S1 subunit which has been proven to be involved 
in receptor binding and host specificity in the well-characterized 
MHV (de Haan et al., 2006). Certain residues in the S protein 
had been associated to enteric or respiratory tropism of BCoV, 
including residues 531, 769 and 1026 (Chouljenko et al., 1998). 
The bubaline strain displayed the same aa residues encountered 
in enteric BCoVs in only one position (D531). The remaining 
two aa residues (S769, G1026) suggest a potential respiratory 
tropism, which is in contrast with the exclusive enteric clini¬ 
cal pattern observed in water buffalo. However, more recent 
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analysis of two enteric and respiratory BCoV pairs failed to 
identify specific aa residue determinants in the S proteins as 
well as in other structural and nonstructural proteins, except for 
a T22I mutation in nsp 4.8 kDa (Zhang et al., 2007). In¬ 
terestingly, strain 179/07-11 showed the same aa determinant 
T22 as the enteric BCoV strains. Other authors have identified a 
truncated form of nsp 4.9 kDa (29 instead of 43 aa) as a genetic 
marker for the respiratory tropism (Gelinas et al., 2001; Vijgen 
et al., 2006). The bubaline strain presented an intact nsp 4.9 kDa 
as enteric BCoVs, but even this aa deletion has been recently 
shown not to be involved in tissue tropism (Zhang et al., 1994). 

Based on the unique biological properties and the more distant 
relatedness of strain 179/07-11 to BCoV with respect to other 
ruminant CoVs, we propose to designate this strain as prototype 
of a host-range variant of BCoV, namely bubaline CoV (BuCoV). 
With regard to the its possible origin, BuCoV may have arisen 
through interspecies transmission of a BCoV strain from cattle to 
water buffaloes, although no recent contacts with cattle or other 
ruminants were reported. This hypothesis is supported by the high 
genetic relatedness between BuCoV and BCoV. Moreover, it 
should be considered that a bovine origin has been strongly 
suggested for other group 2 CoVs less genetically related to 
BCoV, such as HCoV-OC43 (Vijgen et al., 2005), PHEV (Vijgen 
et al., 2006) and CRCoV (Erles et al., 2007). 

Accumulation of point mutations, as well as small insertions 
and deletions in coding and non-coding sequences, are the 
dominant forces in the microevolution of plus-stranded (+) RNA 
viruses, resulting in proliferation of virus strains, serotypes and 
subtypes (Dolja and Carrington, 1992). Coronaviruses are thought 
to mutate at high frequency like most RNA viruses as a con¬ 
sequence of high error rates of the RNA polymerase that are 
predicted to accumulate several base substitutions per round of 
replication (Jarvis and Kirkegaard, 1991; Lai and Holmes, 2001). 
Changes in tissue tropisms and/or interspecies transmission of 
CoVs occur through genetic variations in structural and/or non¬ 
structural proteins (Laude et al., 1993; Vennema et al., 1998; Guan 
et al., 2003; Rottier et al., 2005; Song et al., 2005; Vijgen et al., 
2005; Decaro et al., 2007c). Genetic determinants that may have 
caused the cattle-to-buffalo interspecies transmission have not 
been identified so far. Although this study has detected a CoV 
strain in buffalo calves with severe diarrhea, the pathogenicity of 
this virus and its etiologic role in enteric disease of water buffalo 
still have to be determined more definitively. In addition, epide¬ 
miological studies could assess whether BuCoV is widespread 
among water buffalo herds as well as whether cross-species 
transmission between buffalo and cattle occurs mainly in areas 
where both closely related ruminant species are raised intensively. 

Materials and methods 

(RT-)PCR assay for screening for bovine viral pathogens 

Nucleic acids for (RT-)PCR assays were purified using the 
commercial kits DNeasy Tissue Kit (QIAGEN S.p.A., Milan, 
Italy) and QIAamp® Viral RNA Mini Kit (Qiagen S.p.A., Milan, 
Italy) from the fecal samples and intestinal contents and QIAamp® 
RNeasy Mini Kit (Qiagen S.p.A.) from the tissue samples. 


RNA extracts were used for detection of BCoV (Erles et al., 
2003), toroviruses (Hoet et al., 2002), rotaviruses (Gouvea et al., 
1994), caliciviruses (Jiang et al., 1999), bovine viral diarrhea 
virus (BVDV) (Sullivan and Akkina, 1995) and bovine 
respiratory syncytial virus (Valarcher et al., 1999). Detection 
of bovine herpesvirus types 1 (BoHV-1) (Vilcek, 1993) and 4 
(BoHV-4) (Boemer et al., 1999) was carried out on the DNA 
templates. RT-PCR and PCR assays were performed using 
Superscript™ One-Step RT-PCR for Long Templates (Life 
Technologies) and LA PCR Kit Ver. 2.1 (TaKaRa Bio Inc., 
Shiga, Japan), respectively. 

Bacteriological and parasitological investigations 

The samples were examined for bacterial and parasitic 
pathogens by standardized methods. For bacteriological in¬ 
vestigations, the fecal samples were plated out on 5% sheep 
blood agar and cultured aerobically at 3 7 °C for 24 h to exclude the 
presence of aerobic pathogens. One hundred milligrams of fecal 
samples was resuspended in 900 pi of fluid thioglycolate me¬ 
dium (FTG). Ten-fold dilutions (10” 2 to 10” 8 ) were subsequently 
plated onto 5% sheep blood agar and egg yolk agar with d- 
cycloserine 400 pg/ml. Bacteria were allowed to grow overnight 
at 37 °C in anaerobic condition. Detection of the most common 
enteric parasites was achieved using zinc sulfate flotation. The 
Ziehl Nielsen staining was also performed on fecal samples or 
intestinal sections for detection of Cryptosporidium spp. 

Real-time RT-PCR for quantification of bovine-like CoVs 

A real-time RT-PCR assay based on TaqMan technology 
(G. Elia et al., manuscript in preparation) was used to quantify 
the viral load in samples tested positive for bovine-like CoVs by 
conventional RT-PCR (Erles et al., 2003). The TaqMan assay 
had been proven to be able to detect BCoV as well as the closely 
related CRCoV (Elia et al., manuscript in preparation). BCoV-like 
RNA copy numbers were calculated on the basis of the standard 
curves generated by 10-fold dilutions of a synthetic RNA obtained 
by in vitro transcription of a plasmid containing the M gene of 
BCoV strain 339/07 (Decaro et al., 2007b). Reverse transcription 
was carried out using GeneAmp® RNA PCR (Applied Biosys¬ 
tems, Applera Italia, Monza, Italy), following the manufacturer’s 
recommendations. The quantitative assay targeting the M gene 
was conducted in a 50-pl reaction mixture containing 25 pi of 
IQ™ Supermix (Bio-Rad Laboratories Sri), 600 nM of primers 
BCoV-F (CCTTCATATCTATACACATCAAGTTGTT) and 
BCoV-R (ACCAGCCATTTTAAATCCTTCA), 200 nM of 
probe BCoV-Pb (6FAM“ CCTTCATATCTATACACAT- 
CAAGTTGTT-BHQ1) and 20 pi of c-DNA. The thermal profile 
consisted of activation of iTaq DNA polymerase at 95 °C for 
10 min followed by 45 cycles of denaturation at 95 °C for 15 s and 
annealing extension at 60 °C for 1 min. 

Virus isolation 

Virus isolation was carried out on the buffalo fecal sample 
containing the highest RNA titer of BCoV-like coronavirus 
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(strain 179/07-11), as determined by real-time RT-PCR. Feces 
were homogenized in Dulbecco’s minimal essential medium 
(D-MEM) containing antibiotics (penicillin 5000 IU/ml, strep¬ 
tomycin 2500 pg/ml, amphotericin B 10 pg/ml). After cen¬ 
trifugation at 3000xg for 15 min, the supernatant was used to 
inoculate confluent monolayers of HRT-18 and MDBK cells in 
the presence of trypsin (5 pg/ml). Viral growth was monitored 
by an immunofluorescence (IF) assay using a BCoV-positive 
bovine serum and a rabbit anti-bovine IgG conjugated with 
fluorescein isothiocyanate (Sigma Aldrich srl, Milan, Italy). 

HA and RDE activity 

Two-fold dilutions (starting from dilution 1:2) of the super¬ 
natant of the CoV-infected HRT-18 cells tested positive by the IF 
assay were made in phosphate buffered saline (PBS, pH 7.2) using 
96-well V-plates. Mouse (0.8%) or chicken (0.4%) erythrocytes 
were added to each dilution. Results were read after 1 h at 4 °C or 
37 °C and expressed as the reciprocal of the highest dilution 
producing HA. The plates were then incubated for additional 4 h at 
37 °C to assess RDE activity, which was expressed as the 
reciprocal of the highest dilution of virus resulting in complete 
disappearance of HA (Hasoksuz et al., 1999). 

PCR amplifications of the 3' end of the BCoV-like virus 

The 3' end of the genome of strain 179/07-11 was amplified 
from the original fecal sample using Superscript™ One-Step 
RT-PCR for Long Templates (Life Technologies, Invitrogen. 


Milan, Italy), according to the manufacturer’s instructions. 
Thirteen partially overlapping fragments encompassing from the 
very 3' end of ORFlb to the 5' end of the untranslated region 
(UTR) were amplified using primer pairs designed on conserved 
regions among BCoV-like group 2 coronaviruses (Table 6). 

Sequence analysis and phytogeny 

The PCR-amplified products were sequenced by Genome 
Express (Meylan, France) and the obtained sequences were 
assembled and analyzed using the BioEdit software package 
(Hall, 1999) and the NCBI’s (htttp://www.ncbi.nlm.nih.gov) 
and EMBL’s (http://www.ebi.ac.uk) analysis tools. Additional 
RT-PCR assays and sequencing attempts were performed to 
close gaps between assembled contigs using strain-specific 
primers. ORFs contained in the amplified genomic region were 
determined either with the ORF Finder tool of NCBI or on the 
basis of the similarity to known coronavirus proteins. The ORFs 
identified in this manner were translated and the predicted 
amino acid (aa) sequences were saved as individual files for 
further analyses. Putative N-glycosylation sites were predicted 
using the NetNglyc 1.0 Server program (http://www.cbs.dtu.dk/ 
services/NetNGlyc/), whereas the NetOGlyc 3.1 Server program 
(http://www.cbs.dtu.dk/services/NetOGlyc/) was used for predic¬ 
tion of O-linked glycosylation sites. Phylogenetic and molecular 
evolutionary analyses were conducted using Mega3 (Kumar et al., 
2004). Phylogenetic trees, based on the nonstmctural 32-kDa 
protein, on the structural spike (S), envelope (E), membrane (M) 
and nucleocapsid (N) proteins and on the genomic RNA 3' end 


Table 6 


Primers used for RT-PCR amplification and sequence analysis 


Primer 

Sequence 5' to 3' 

Sense 

Position a 

Amplicon size (bp) 

BCV-22001F 

TAGACTTGAAATAGTTAAGCTTGGTG 

+ 

22,001-22,026 

893 

BCV-22894R 

A A ATTAGCTTC ACGAGCTATATAT GC 

— 

22,868-22,893 


BCV-22769F 

TATCGCAGCCTTACTTTTGTTAATG 

+ 

22,769-22,793 

527 

BCV-23295R 

CGAAAATAACAGTACGGGGGTTGACA 

— 

23,270-23,295 


BCV-23112F 

TACCCTCTGGTAATTATTTAGCCATTTCA 

+ 

23,112-23,140 

776 

BCV-23887R 

TTCCCTTCAGTGCCATATTACGATATGT 

— 

23,860-23,887 


BCV-23510F 

TAT GAT C CGCTAC C A ATTATTTT GCTT GGC A 

+ 

23,510-23,540 

817 

BCV-24326R 

AC A AC ACC AGT GTCTGTA A AATAT GCA 

— 

24,300-24,326 


BCV-24182F 

TAG A ACTAT GGC ATT GGG ATAC AGGT GTT G 

+ 

24,182-24,211 

1254 

BCV-25435R 

TACACCTATCCCCTTGTAAACAAGAGTA 

— 

25,409-25,435 


BCV-25301F 

ACTTAGTT GGC ATAGGT G AGC ACTGTT C 

+ 

25,301-25,328 

851 

BCV-26151R 

ACATGCTACATAATCACCACAGACAA 

— 

26,126-26,151 


BCV-26028F 

TTT GTAT G A A ATT C A A ATAC CTT C AG AG 

+ 

26,028-26,055 

877 

BCV-26904R 

GTCTATCTGAGCTTGCGCTTCAAGAGCA 

— 

26,877-26,904 


BCV-26760F 

TA A A ATTC A AGCT GTT GTTA AT GCA A AT 

+ 

26,760-26,787 

1059 

BCV-27818R 

GCTCGACCTAAATGGGTCTTATAATTAGA 

— 

27,791-27,818 


BCV-27667F 

GGT GGTT GTT GT GAT G ATTATAC T GG AC A 

+ 

27,667-27,695 

853 

BCV-28519R 

ACTACAACTATTATAACCAATAAACACAT 

— 

28,491-28,519 


BCV-28380F 

TCTTAGCTGTTGACTTTATTACCTGG 

+ 

28,380-28,405 

919 

BCV-29298R 

ACATAAACAGCAAAACCACTAGTATCGCC 

— 

29,270-29,298 


BCV-29160F 

AAGGTATAAAACTAGGTACTGGCTAT 

+ 

29,160-29,185 

876 

BCV-30035R 

TGCGCGATCCTGCACTAGAGGCTCTAC 

— 

30,009-30,035 


BCV-29900F 

GAGGCTATTCCGACTAGGTTTCCGCCT 

+ 

29,900-29,926 

827 

BCV-30726R 

GTGTCTTCAGTAAAGGGCTCATCCATC 

— 

30,700-30,726 


BCV-30610F 

TGATAATATAAGTGTTGCAGCGCCCAAA 

+ 

30,610-30,637 

419 

BCV-31028R 

GTGATTCTTCCAATTGGCCATAATT 

— 

31,004-31,028 



a Primers position is referred to the sequence of BCoV strain ENT (accession number AF391541). 
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(9.6 kb) of CoV strain 179/07-11, were elaborated using both 
parsimony and neighbor-joining methods, supplying a statistical 
support with bootstrapping over 1000 replicates. 

SimPlot analysis 

The SimPlot program (version 3.2) was used to analyze the 
genetic distance of the 3' end of the genomes of ruminant Co Vs, 
including the bubaline strain, in reference to the same region of 
BCoV strain ENT and this genetic distance was plotted versus 
nucleotide (nt) positions (Lole et al., 1999). 

Nucleotide sequence accession number 

The nt sequence of the bubaline strain 179/07-11 has been 
deposited in GenBank under accession number EU019216. 

Acknowledgments 

This work was supported by grants from University of Bari, 
Italy: project ex 60% 2007 “Messa a punto di un sistema real-time 
RT-PCR per la identificazione e la quantificazione dell’RNA del 
coronavirus bovino”. 

References 

Benfield, D.A., Saif, L.J., 1990. Cell culture propagation of a coronavirus 
isolated from cows with winter dysentery. J. Clin. Microbiol. 28,1454-1457. 
Boemer, B., Weigelt, W., Buhk, H.J., Castrucci, G., Ludwig, H., 1999. A sensitive 
and specific PCR/Southem blot assay for detection of bovine herpesvirus 4 in 
calves infected experimentally. J. Virol. Methods 83, 169-180. 

Cho, K.O., Halbur, P.G., Bmna, J.D., Sorden, S.D., Yoon, K.J., Janke, B.H., 
Chang, K.O., Saif, L.J., 2000. Detection and isolation of coronavirus from 
feces of three herds of feedlot cattle during outbreaks of winter dysentery- 
like disease. J. Am. Vet. Med. Assoc. 217, 1191-1194. 

Chouljenko, V.N., Kousoulas, K.G., Lin, X., Storz, J., 1998. Nucleotide and 
predicted amino acid sequences of all genes encoded by the 3' genomic 
portion (9.5 kb) of respiratory bovine coronaviruses and comparisons among 
respiratory and enteric coronavimses. Virus Genes 17, 3-42. 

Chouljenko, V.N., Lin, X.Q., Storz, J., Kousoulas, K.G., Gorbalenya, A.E., 
2001. Comparison of genomic and predicted amino acid sequences of 
respiratory and enteric bovine coronaviruses isolated from the same animal 
with fatal shipping pneumonia. J. Gen. Virol. 82, 2927-2933. 
de Haan, C.A., Te Lintelo, E., Li, Z., Raaben, M., Wurdinger, T., Bosch, B.J., 
Rottier, P.J., 2006. Cooperative involvement of the SI and S2 subunits of the 
murine coronavirus spike protein in receptor binding and extended host 
range. J. Virol. 80, 10909-10918. 

Decaro, N., Desario, C., Elia, G., Mari, V., Lucente, M.S., Cordioli, P., 
Colaianni, M.L., Martella, V., Buonavoglia, C., 2007a. Serological and 
molecular evidence that canine respiratory coronavirus is circulating in Italy. 
Vet. Microbiol. 121, 225-230. 

Decaro, N., Mari, V., Desario, C., Campolo, M., Elia, G., Martella, V., Greco, G., 
Cirone, F., Colaianni, M.L., Cordioli, P., Buonavoglia, C., 2007b. Severe 
outbreak of bovine coronavirus infection in dairy cattle during the warmer 
season. Vet. Microbiol., doi:10.1016/j.vetmic.2007.06.024 Electronic publi¬ 
cation ahead of print. 

Decaro, N., Martella, V., Elia, G., Campolo, M., Desario, C., Cirone, F., 
Tempesta, M., Buonavoglia, C., 2007c. Molecular characterisation of the 
virulent canine coronavirus CB/05 strain. Virus Res. 125, 4-60. 

Decaro, N., Campolo, M., Desario, C., Cirone, F., D’abramo, M., Lorusso, E., 
Greco, G., Mari, V., Colaianni, M.L., Elia, G., Martella, V., Buonavoglia C., 
in press. Respiratory disease associated with bovine coronavirus infection in 
cattle herds in Southern Italy. J. Vet. Diagn. Invest. 


Dolja, V.V., Carrington, J.C., 1992. Evolution of positive-strand RNA viruses. 
Semin. Virol. 3, 315-326. 

Enjuanes, L., Brian, D., Cavanagh, D., Holmes, K., Lai, M.M.C., Laude, H., 
Masters, P., Rottier, P., Siddell, S., Spaan, W.J.M., Taguchi, F., Talbot, R, 
2000. Family Coronaviridae. In: van Regenmortel, M.H.V., Fauquet, C.M., 
Bishop, D.H.L., Carstens, E.B., Estes, M.K., Lemon, S.M., Maniloff, J., 
Mayo, M.A., McGeoch, D.J., Pringle, C.R., Wickner, R.B. (Eds.), Virus 
Taxonomy, Classification and Nomenclature of Viruses. Academic Press, 
New York, pp. 835-849. 

Erles, K., Toomey, C., Brooks, H.W., Brownlie, J., 2003. Detection of a group 2 
coronavirus in dogs with canine infectious respiratory disease. Virology 310, 
216-223. 

Erles, K., Shiu, K.B., Brownlie, J., 2007. Isolation and sequence analysis of 
canine respiratory coronavirus. Virus Res. 124, 78-87. 

Fukutomi, T., Tsunemitsu, H., Akashi, H., 1999. Detection of bovine 
coronaviruses from adult cows with epizootic diarrhea and their antigenic 
and biological diversities. Arch. Virol. 144, 997-1006. 

Gelinas, A.M., Boutin, M., Sasseville, A.M., Dea, S., 2001. Bovine coronavimses 
associated with enteric and respiratory diseases in Canadian dairy cattle display 
different reactivities to anti-HE monoclonal antibodies and distinct amino acid 
changes in their HE, S and ns4.9 protein. Vims Res. 76, 43-57. 

Gouvea, V., Santos, N., Mdo, C., 1994. Identification of bovine and porcine 
rotavims G types by PCR. J. Clin. Microbiol. 32, 1338-1340. 

Guan, Y., Zheng, B.J., He, Y.Q., Liu, X.L., Zhuang, Z.X., Cheung, C.L., Luo, S.W., 
Li, P.H., Zhang, L.J., Guan, Y.J., Butt, K.M., Wong, K.L., Chan, K.W., Lim, 
W., Shortridge, K.F., Yuen, K.Y., Peiris, J.S., Poon, L.L., 2003. Isolation and 
characterization of viruses related to the SARS coronavims from animals in 
southern China. Science 302, 276-278. 

Guy, J.S., Breslin, J.J., Breuhaus, B., Vivrette, S., Smith, L.G., 2000. 
Characterization of a coronavirus isolated from a diarrheic foal. J. Clin. 
Microbiol. 38, 4523-4526. 

Hall, T.A., 1999. BioEdit: a user-friendly biological sequence alignment and 
analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41,95-98. 

Hasoksuz, M., Lathrop, S.L., Gadfield, K.L., Saif, L.J., 1999. Isolation of 
bovine respiratory coronavimses from feedlot cattle and comparison of their 
biological and antigenic properties with bovine enteric coronavimses. Am. J. 
Vet. Res. 60, 1227-1233. 

Hasoksuz, M., Alekseev, K., Vlasova, A., Zhang, X., Spiro, D., Halpin, R., Wang, 
S., Ghedin, E., Saif, L.J., 2007. Biologic, antigenic, and full-length genomic 
characterization of a bovine-like coronavims isolated from a giraffe. J. Virol. 
81,4981-4990. 

Hoet, A.E., Cho, K.O., Chang, K.O., Loerch, S.C., Wittum, T.E., Saif, L.J., 
2002. Enteric and nasal shedding of bovine torovims (Breda vims) in feedlot 
cattle. Am. J. Vet. Res. 63, 342-348. 

Jarvis, T.C., Kirkegaard, K., 1991. The polymerase in its labyrinth: mechanisms 
and implications of RNA recombination. Trends Genet. 7, 186-191. 

Jiang, X., Huang, P.W., Zhong, W.M., Farkas, T., Cubitt, D.W., Matson, D.O., 
1999. Design and evaluation of a primer pair that detects both Norwalk- and 
Sapporo-like calicivimses by RT-PCR. J. Virol. Methods 83, 145-154. 

Jin, L., Cebra, C.K., Baker, R.J., Mattson, D.E., Cohen, S.A., Alvarado, D.E., 
Rohrmann, F., 2007. Analysis of the genome sequence of an alpaca 
coronavims. Virology 365, 198-203. 

Kumar, S., Tamura, K., Nei, M., 2004. MEGA3: integrated software for 
Molecular Evolutionary Genetics Analysis and sequence alignment. Brief. 
Bioinform. 5, 150-163. 

Lai, M.M.C., Holmes, K.V., 2001. Coronaviridae: the vimses and their 
replication, In: Knipe, D.M., Howley, P.M. (Eds.), Fields Virology, 4th 
edition. Lippincott Williams and Wilkins, Philadelphia, PA, pp. 1163-1185. 

Lathrop, S.L., Wittum, T.E., Brock, K.V., Loerch, S.C., Perino, L.J., Bingham, 
H.R., McCollum, F.T., Saif, L.J., 2000. Association between infection of the 
respiratory tract attributable to bovine coronavims and health and growth 
performance of cattle in feedlots. Am. J. Vet. Res. 61, 1062-1066. 

Laude, H., Van Reeth, K., Pensaert, M., 1993. Porcine respiratory coronavims: 
molecular features and vims-host interactions. Vet. Res. 24, 125-150. 

Lole, K.S., Bollinger, R.C., Paranjape, R., Gadkari, S., Kulkami, D., Novak, S.S., 
Ingersoll, N.G., Sheppard, R., Ray, H.W., 1999. Full-length human im¬ 
munodeficiency vims type 1 genomes from subtype C-infected seroconverters 
in India, with evidence of intersubtype recombination. J. Virol. 73, 152-160. 


222 


N. Decaro et al. / Virology 370 (2008) 213-222 


Muniiappa, L., Mitov, B.K., Kharalambiev, Kh.E., 1985. Demonstration of 
coronavirus infection in buffaloes. Vet. Med. Nauki. 22, 27-32. 

Romano, R, Ricciardi, A., Salzano, G., Suzzi, G., 2001. Yeasts from water 
buffalo mozzarella, a traditional cheese of the Mediterranean area. Int. J. 
Food Microbiol. 69, 45-51. 

Rottier, P.J., Nakamura, K., Schellen, R, Volders, H., Haijema, B.J., 2005. 
Acquisition of macrophage tropism during the pathogenesis of feline 
infectious peritonitis is determined by mutations in the feline coronavirus 
spike protein. J. Virol. 79, 14122-14130. 

Saif, L.J., Brock, K.V., Redman, D.R., Kohler, E.M., 1991. Winter dysentery in 
dairy herds: electron microscopic and serological evidence for an association 
with coronavirus infection. Vet. Rec. 128, 447-449. 

Saif, L.J., Redman, D.R., Brock, K.V., Kohler, E.M., Heckert, R.A., 1998. 
Winter dysentery in adult dairy cattle: detection of coronavirus in the faeces. 
Vet. Rec. 123, 300-301. 

Snodgrass, D.R., Terzolo, H.R., Sherwood, D., Campbell, I., Menzies, J.D., Synge, 
B.A., 1986. Aetiology of diarrhoea in young calves. Vet. Rec. 119, 31-34. 

Song, H.D., Tu, C.C., Zhang, G.W., Wang, S.Y., Zheng, K., Lei, L.C., Chen, Q.X., 
Gao, Y.W., Zhou, H.Q., Xiang, H., Zheng, H.J., Chern, S.W., Cheng, F., 
Pan, C.M., Xuan, Eh, Chen, S.J., Luo, H.M., Zhou, D.H., Liu, Y.F., He, J.F., 
Qin, P.Z., Li, L.H., Ren, Y.Q., Liang, W.J., Yu, Y.D., Anderson, L., Wang, 
M., Xu, R.H., Wu, X.W., Zheng, H.Y., Chen, J.D., Liang, G., Gao, Y., Liao, 
M., Fang, L., Jiang, L.Y., Li, H., Chen, F., Di, B., He, L.J., Lin, J.Y., Tong, 
S., Kong, X., Du, L., Hao, P., Tang, H., Bernini, A., Yu, H.J., Spiga, O., 
Guo, Z.M., Pan, H.Y., He, W.Z., Manuguerra, J.C., Fontanet, A., Danchin, 
A., Niccolai, N., Li, Y.X., Wu, C.I., Zhao, G.P., 2005. Cross-host evolution 
of severe acute respiratory syndrome coronavirus in palm civet and human. 
Proc. Natl. Acad. Sci. U. S. A. 102, 2430-2435. 

Storz, J., Purdy, C.W., Lin, X., Burrell, M., Truax, R.E., Briggs, R.E., Frank, G.H., 
Loan, R.W., 2000. Isolation of respiratory bovine coronavirus, other cytocidal 
vimses, and Pasteurella spp from cattle involved in two natural outbreaks of 
shipping fever. J. Am. Vet. Med. Assoc. 216, 1599-1604. 

Sullivan, D.G., Akkina, R.K., 1995. A nested polymerase chain reaction assay to 
differentiate pestiviruses. Virus Res. 38, 231-239. 

Traven, M., Naslund, K., Linde, N., Linde, B., Silvan, A., Fossum, C., Hedlund, 
K.O., Larsson, B., 2001. Experimental reproduction of winter dysentery in 


lactating cows using BCV—comparison with BCV infection in milk-fed 
calves. Vet. Microbiol. 81, 127-151. 

Tsunemitsu, H., el-Kanawati, Z.R., Smith, D.R., Reed, H.H., Saif, L.J., 1995. 
Isolation of coronaviruses antigenically indistinguishable from bovine 
coronavirus from wild ruminants with diarrhea. J. Clin. Microbiol. 33, 
3264-3269. 

Valarcher, J.F., Bourhy, H., Gelfi, J., Schelcher, F., 1999. Evaluation of a nested 
reverse transcription-PCR assay based on the nucleoprotein gene for diagnosis 
of spontaneous and experimental bovine respiratory syncytial vims infections. 
J. Clin. Microbiol. 37, 1858-1862. 

Vennema, H., Poland, A., Foley, J., Pedersen, N.C., 1998. Feline infectious 
peritonitis vimses arise by mutation from endemic feline enteric corona- 
viruses. Virology 243, 150-157. 

Vijgen, L., Keyaerts, E., Moes, E., Thoelen, I., Wollants, E., Lemey, P., 
Vandamme, A.M., Van Ranst, M., 2005. Complete genomic sequence of 
human coronavims OC43: molecular clock analysis suggests a relatively 
recent zoonotic coronavims transmission event. J. Virol. 79, 1595-1604. 

Vijgen, L., Keyaerts, E., Lemey, P., Maes, P., Van Reeth, K., Nauwynck, H., 
Pensaert, M., Van Ranst, M., 2006. Evolutionary history of the closely 
related group 2 coronavimses: porcine hemagglutinating encephalomyelitis 
vims, bovine coronavims, and human coronavims OC43. J. Virol. 80, 
7270-7274. 

Vilcek, S., 1993. Detection of the bovine herpesvirus-1 (BHV-1) genome by 
PCR. J. Virol. Methods 41, 245-247. 

Woo, P.C., Lau, S.K., Chu, C.M., Chan, K.H., Tsoi, H.W., Huang, Y., Wong, B.H, 
Poon, R.W., Cai, J.J., Luk, W.K., Poon, L.L., Wong, S.S., Guan, Y., Peiris, J.S., 
Yuen, K.Y., 2005. Characterization and complete genome sequence of a novel 
coronavims, coronavims HKU1, from patients with pneumonia. J. Virol. 79, 
884-895. 

Zhang, X.M., Herbst, W., Kousoulas, K.G., Storz, J., 1994. Biological and 
genetic characterization of a hemagglutinating coronavims isolated from a 
diarrhoeic child. J. Med. Virol. 44, 152-161. 

Zhang, X., Hasoksuz, M., Spiro, D., Halpin, R., Wang, S., Vlasova, A., Janies, 
D., Jones, L.R., Ghedin, E., Saif, L.J., 2007. Quasispecies of bovine enteric 
and respiratory coronavimses based on complete genome sequences and 
genetic changes after tissue culture adaptation. Virology 363, 1-10. 


